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Abstract 

This paper presents an analysis of the inclusive properties of diffractive deep inelastic 
scattering events produced in ep interactions at HERA. The events are characterised by 
a rapidity gap between the outgoing proton system and the remaining hadronic system. 
Inclusive distributions are presented and compared with Monte Carlo models for diffractive 
processes. The data are consistent with models where the pomeron structure function has 
a hard and a soft contribution. The diffractive structure function is measured as a function 
of Xjp, the momentum fraction lost by the proton, of /3, the momentum fraction of the 
struck quark with respect to Xp, and of Q^. The Xp dependence is consistent with 
the form where a = 1.30 ± 0.08 {stat) t oaI (sys) in ah bins of (3 and Q^. 

In the measured range, the diffractive structure function approximately scales with 
at fixed (3. In an Ingelman-Schlein type model, where commonly used pomeron flux 
factor normalisations are assumed, it is found that the quarks within the pomeron do not 
saturate the momentum sum rule. 
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1 Introduction 



We present an analysis of deep inelastic scattering (DIS) events with a large rapidity gap 
between the outgoing proton system and the remaining hadronic final state. The general 
properties of these events indicate that the underlying production mechanism is leading twist 
and diffractive ^ Q]. Diffractive processes are generally understood to proceed through the 
exchange of a colourless object with the quantum numbers of the vacuum, generically called 
the pomeron [0]. The true nature of this exchanged "object" remains unclear. 

The analysis of soft hadron-hadron collisions implies that pomeron-exchange can be described 
by a pomeron-hadron coupling constant and a pomeron-propagator |Q . This led to the propo- 
sition of Ingelman and Schlein ^ to treat the pomeron as a quasi-real particle which is emitted 
by a hadron, described in terms of a parton density and characterised by a structure function 
F,^ which can be studied in deep inelastic scattering. The assumption of factorisation implies 
that the pomeron structure is independent of the process of emission. Evidence for a partonic 
structure of the pomeron was observed by the UA8 collaboration ^ and later by the HERA 
experiments p, ^, ||, ^ . These data also gave a first insight into the structure function of the 
pomeron. The UA8 data show a predominantly hard structure, where on average the par- 
tons carry a large fraction of the momentum of the pomeron. However, these data could not 
distinguish between the quark and the gluon content of the pomeron. 

This paper presents a study of the structure of the pomeron in DIS at HERA. We discuss 
first the variables and cuts which are used to isolate diffractive DIS events. The observed 
kinematic distributions of the data are compared to Monte Carlo models of diffractive processes 
in a region where the diffractive contribution dominates. A measurement of the diffractive 
structure function is presented, integrated over t, the square of the momentum transfer at 
the proton vertex, as a function of Xp, the momentum fraction lost by the proton, of /3, 
the momentum fraction of the struck quark with respect to x^, and of Q^. The data are 
used to determine the dependence at fixed (3 in order to test factorisation; extract the P 
dependence of the diffractive structure function at fixed Q^; investigate the dependence of 
the diffractive structure function at fixed (3, in order to test scale invariance; and, examine the 
general dependence onx^, (3 and by comparing the data with different models for diffractive 
dissociation. 



2 Experimental setup 
2.1 HERA 

This analysis is based on data collected with the ZEUS detector at the electron-proton collider 
HERA. During 1993, HERA was operated at a proton energy, Ep, of 820 GeV and an electron 
energy, Ee, of 26.7 GeV. HERA is designed to run with 210 bunches in each of the electron and 
proton rings, with an interbunch spacing of 96 ns. For the 1993 data-taking 84 paired bunches 
were filled for each beam and in addition 10 electron and 6 proton bunches were left unpaired 
for background studies. Typical total currents were 10 mA for both beams. 
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2.2 The ZEUS detector 



Details of the ZEUS detector can be found in |TD, [T^. The following is hence restricted to a 
short description of the components relevant to the present analysis. 

Charged particles are tracked by the inner tracking detectors which operate in a magnetic field 
of 1.43 T provided by a thin superconducting coil. Immediately surrounding the beampipe is 
the vertex detector (VXD) which consists of 120 radial cells, each with 12 sense wires ||12|| . 
Surrounding the VXD is the cylindrical central tracking detector (CTD) which consists of 72 
cylindrical drift chamber layers, organised into 9 superlayers |jl3|] (5 axial and 4 small angle 
stereo layers). In events with charged tracks, using the combined data from both chambers, 
resolutions of 0.4 cm in Z and 0.1 cm in radius in the XY planeQ are obtained for the primary 
vertex reconstruction. 

The high resolution uranium-scintillator calorimeter (CAL) |T^, |T3|, ^ consists of three parts, 
forward (FCAL) covering the pseudorapidityf] region 4.3 > ?7 > 1.1, barrel (BCAL) covering 
the central region 1.1 > r/ > —0.75 and rear (RCAL) covering the backward region —0.75 > 
ri > —3.8. Holes of 20 x 20 cm^ in the centre of FCAL and RCAL are required to accommodate 
the HERA beam pipe. The resulting solid angle coverage is 99.7% of An. The calorimeter 
is subdivided longitudinally into one electromagnetic (EMC) section and one (RCAL) or two 
(FCAL and BCAL) hadronic (HAC) sections. The sections are subdivided into cells, each of 
which is read out by two photomultiplier tubes. The CAL also provides a time resolution 
of better than 1 ns for energy deposits greater than 4.5 GeV, which is used for background 
rejection. 

The C5 beam monitor, a small lead-scintillator counter assembly around the beam pipe located 
at Z = —3.2 m, has been used to measure the timing and longitudinal structure of the proton 
and electron bunches, and to reject events from upstream proton-gas interactions. The vetowall 
detector, consisting of two layers of orthogonal scintillator strips on either side of an 87 cm thick 
iron wall centred aX Z = —7.3 m, was also used to tag upstream background events. 

The luminosity is measured from the rate observed in the luminosity photon detector of hard 
bremsstrahlung photons from the Bethe-Heitler process ep —>■ e'pj. The luminosity detector 



ITf consists of a photon and an electron lead-scintillator calorimeter. Bremsstrahlung photons 
emerging from the electron-proton interaction point at angles below 0.5 mrad with respect to 
the electron beam axis hit the photon calorimeter placed at 107 m distance along the electron 
beam line. Electrons emitted at scattering angles less than 5 mrad and with energies 0.2Ee < 
E'^ < 0.9Ee are deflected by beam magnets and hit the electron calorimeter placed 35 m from 
the interaction point. 



2.3 Trigger conditions 

Data were collected with a three level trigger |TD|; details of the first (FLT) and second (SLT) 



level decision for DIS events can be found in previous publications |jT 

^The ZEUS coordinate system is defined as right handed with the Z axis pointing in the proton beam 
direction, hereafter referred to as forward, the X axis pointing towards the centre of HERA and the Y axis 
pointing upwards. 

^The pseudorapidity r] is defined as — ln(tan |), where the polar angle 9 is taken with respect to the proton 
beam direction and, in this case, refers to the nominal interaction point. 
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The third level trigger (TLT) rejects beam-gas events using timing and cosmic-ray events by 
a combination of timing and topology. Finally it passes the accepted events through a set of 
filters in order to categorise the events. The category of DIS neutral current events is defined 
by the requirement of an electron candidate in the RCAL or BCAL. A cut was performed on 
6 = Eji?j(l — cosOi) > 20GeV — 2E^, where Ei and 6i are the energies and polar angles (in this 
case, with respect to the nominal interaction point) of calorimeter cells and is the energy 
measured in the photon calorimeter of the luminosity monitor. For events fully contained in the 
main detector 6 ~ 2Ee = 53.4 GeV, whereas for low-Q^ events the scattered electron escapes 
through the rear beam pipe and S peaks at low values. 



3 Kinematic variables 

The kinematic variables used to describe DIS events 

e (k) + p (P) — > e' {k') + anything 

are the following: the negative of the squared four-momentum transfer carried by the virtual 
photonQ: 

Q' = -q' = -{k - k'f- 

the Bjorken variable: 

X 



2P-q 

the variable which describes the energy transfer to the hadronic final state: 

P ■ q 

y=T-k' 

and the centre-of-mass energy W of the virtual-photon proton (7*p) system, where: 

= {q + Pf = ^^^^^—^ + Ml 

X 

with Mp denoting the proton mass. 

These variables, only two of which are independent at fixed ep centre-of-mass energy squared 
s = {k + P)'^, can be reconstructed in a variety of ways using combinations of electron and 
hadronic system energies and angles [|19|. The variable y, calculated from the electron variables, 
is given by: 

^ El 1-cos^; 

Ve = 1 

^ Ee 2 

where E'^, 6'^ denote the energy and angle of the scattered electron. Alternatively, y can be 
estimated from the hadronic system, using the Jacquet-Blondel technique pO| : 



Y^iEiil- cos9i 

yjB = - 



2 ■ K 



•^In the range used for this analysis, ep interactions are described to sufficient accuracy by the exchange 
of a virtual photon. 
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where Ei and 6i are the energies and polar angles of calorimeter cells which are associated with 
the hadronic system. 

Studies of the kinematic variables have shown that for this analysis it is advantageous to use 
the so-called double angle (DA) method, in which the angles of the scattered electron and the 
hadron system are used to determine x and Q^. Quantities determined in this way will be 
denoted by the subscript DA. Formulae to calculate Qjjj^, Wda, xda and i/da are given in |T^ 



In the diffractive DIS process shown in Fig. 1: 

e{k) + p (P) ^ e' (k') + p\P') + X, 

the hadronic system X (exclusive of the proton) and the scattered electron e' are detected in 
the main detector. The proton remnant p' remains undetected. When the system X is fully 
contained its invariant mass, Mx, can be determined from the calorimeter cell information as 
follows Denoting the energy, momentum and polar angle of the final hadronic system as 
Eh, Ph and Oh-, respectively; and pi as the vector constructed from the energy polar angle 
9i and azimuthal angle (pi of cell i\ then: 

cosOh = I ^ (1) 
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Ph- . 2n 

sm Uh 
Eh = Ph cos 9h + 2E^yDA 



from which Mx is determined by the definition Mx = y Ejj — pfj. 
The squared four-momentum transfer at the proton vertex is given by: 

t = {P- P'f, 

whose absolute magnitude is expected to be small compared to Q"^ + Mx in diffractive processes 
for the kinematic region studied here. To describe diffractive deep inelastic scattering, in 
addition to x and Q^, the following variables are used: 

_{P-P')-q _ Ml + Q'^-t ^ Ml + Q"^ 

X r, 



'IP 



/3 



P-q W^ + Q^-M^ + g2 

Q' X Q2 



2{P-P')-q x^ Ml + Q^-t Ml + Q'^' 

In models where diffraction is described by the exchange of a particle-like pomeron, x^ is the 
momentum fraction of the pomeron in the proton and j3 is the momentum fraction of the struck 
quark within the pomeron. For the structure of the pomeron in DIS, the variable (3 plays a role 
analogous to that of Bjorken-x for the structure of the proton. 
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4 DifFractive structure function 



For unpolarised beams, the differential cross section for single diffractive dissociation can be 
described in terms of the diffractive structure function, {f3,Q'^,Xjp,t): 

where a is the electromagnetic coupling constant and the 6i denote corrections due to Z° 
exchange and due to radiative corrections which are small in the measured range. The contri- 
bution of Fl to the diffractive cross section is not known. If such a term were included, the F2 
values would become larger at large y values (corresponding to small values). The effect of 
this uncertainty is considered in section |8?T| . Note that the function F^^'^\j3,Q'^,Xjp,t) can be 
related to that of J-'2^^\x, Q^, Xjp,t). Integrating JF^*-^-* over x^ and t one can directly compare 
it to the inclusive proton structure function F2(x, 



In this analysis, an integral is performed over t, corresponding to the (undetected) momentum 
transfer to the proton system. For this initial measurement we neglect the effect of Fl and 
the additional contributions noted above, yielding the following expression for F^^^^ , where the 
cross section is evaluated as a function of j3, and : 



dfSdQ^dx^ f3Q^ 



[l + {l-yr)Fr'{f3,Q' 



following the procedure of |22|, where the relation x = f5Xjp has been used. 



5 DifFractive models and Monte Carlo simulation 



Different approaches exist to model diffractive processes such as that depicted in Fig. |I[ In 
this paper we compare the data with the predictions of the factorisable models of Ingelman 
and Schlein [Q, Donnachie and Landshoff and Capella et al. ||2^, as well as the non- 

We have earlier found that the Monte 



[Q, Donnachie and Landshoff J 
factorisable model of Nikolaev and Zakharov 
Carlo implementations p 



28[] of the models described above provide reasonable descriptions of 



the shape of the energy flow and of the observed fraction of events with one or more jets p|, ^ 



In the model of Ingelman and Schlein the proton emits a pomeron which is treated as a 
(virtual) hadron whose structure is probed by the virtual photon. The pomeron is described by 
a structure function F^{f3,Q'^) which is independent of the process of emission. In this sense 
factorisation is predicted in the model: 



(AQ' 



X J, 



t) = fip{x^,t)-Ff{P,Q' 



The flux factor f{Xjp,t), describing the flux of pomerons in the proton, can be extracted from 
hadron-hadron scattering with an accuracy of approximately 30%, assuming universality of the 



pomeron flux. A comparison of different flux factors can be found in ||29| 
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For this analysis we used the POMPYT Monte Carlo implementation of the Ingelman- 
Schlein model. Two samples of events were generated, corresponding to a hard quarkonic 
structure function, 

and to a soft quarkonic structure function. 

The two samples are denoted by "Hard Pomeron" (HP) and "Soft Pomeron" (SP) respectively. 
The normalisation constant 5/3 is based on the assumption that the momentum sum rule 
(MSR) is satisfied for two light quark flavours (u,d). If s quarks would have to be included the 
normalisation factor would be reduced from 5/3 to 4/3 ||31[. The dependence is expected 



to be weak and is neglected. The Ingelman-Schlein form of the flux is parametrised by a fit to 
UA4 data ^ 0]: 

fip{x^,t) = l^^^ ■ (6.38 + 0.424 e^^*). 

In the Donnachie-Landshoff (DL) model diffraction in DIS is described through pomeron ex- 
change between the virtual photon and the proton, with the pomeron coupling predominantly 
to quarks |^ . The authors calculate the cross section in the framework of Regge theory. The 



result can be interpreted in terms of a pomeron structure function with the resulting /5 de- 
pendence similar to HP but with a normalisation which is calculated to be approximately a 
factor of 6.2 smaller. The authors also predict an additional soft contribution to the pomeron 
structure function which is expected to become important only for (3 < 0.1. The flux factor. 
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/^(x^,t) = ^Fi(t)^a:/-^"W, 

is related to the elastic form factor of the proton, Fi{t) = "^'^2"^f* ( i_f;/o7 )^; and to the pomeron- 
quark coupling, /?o — 1-8 GeV~^, extracted from hadron-hadron data. The term represents 
the pomeron propagator with the pomeron trajectory, a{t) = 1.085 + 0.25 ■ t. Therefore the 
dependence of F^^'^^ is controlled by the pomeron trajectory (^2^*-''^ oc x^p^"^"'^*^). Integrated 
over t, the predicted effective x^, -dependence at flxed /3 is approximately (l/x^)^'°^ in the 
measured range of x^. 



In a recent publication ||53[ Goulianos proposed to use a modifled flux factor, which is renor- 
malised to unity for flxed centre-of-mass energy W. Also a modiflcation of the t-dependence 



and the pomeron trajectory according to recent CDF data |3J] is proposed: 



iVo 



where a{t) = 1.115 + 0.26 ■ t and A^o is a normalisation factor which can be approximated by 

-400 



w^^o.23 Integrated over t, the effective x^, -dependence of this flux factor is approximately 



[l/x 



,0.93 

p J 



Capella et al. calculate the diffractive structure function in the framework of conventional 
Regge theory p4[. Using Regge factorisation, they relate the pomeron structure function to 
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the deuteron structure function using parameters which are determined from soft hadronic 
diffraction data with an appropriate change for the disappearance of screening corrections with 
increasing Q"^. For Q"^ = 10 GeV^ they obtain: 

Ff Q') = a- . (1 - + 0.015 ■ p-'-^' ■ (1 - , 

where a is estimated to be in the range 0.04 to 0.06. We chose a = 0.06 for comparison with 
the data. 

In the model of Nikolaev and Zakharov diffractive dissociation is described as a fluctuation 



of the photon into a qq or qqg Fock state The interaction with the proton proceeds 



via the exchange of a BFKL type pomeron, starting in lowest-order from the exchange of 
a Low-Nussinov 0] pomeron which corresponds to two gluons in a colour-singlet state. The 
result for the cross section can be approximated by a two-component structure function of the 
pomeron, each component having its own flux factor. This corresponds to factorisation breaking 
which is caused by BFKL evolution effects. The result for the "hard" component reflects the 
case where the photon fluctuates into a qq pair and leads to a (3 dependence similar to that of 
the HP and the DL models with a normalisation closer to the latter (but with very different 
predictions for the contribution of heavy flavours). The "soft" contribution, which reflects the 
case where the photon fluctuates into qqg, is assumed to be proportional to (1 — /?)^ and the 
normalisation is fixed by the triple pomeron coupling. The relative size of these contributions 
and the overall normalisation are predicted with an uncertainty of about 30%. evolution 
effects have been calculated for this model, and have been found to be rather weak for (3 > 0.1. 
We used a Monte Carlo implementation of this model [EH] which is based on the cross section 



given in and is interfaced to the Lund fragmentation scheme. We refer to this model as 



NZ. In this implementation the mass spectrum contains both components but the qqg states 
are fragmented into hadrons as if they were a qq system with the same Mx. 

Both Monte Carlo generators have limitations in the generation of small masses Mx 
{Mx < 1.7 GeV for NZ and Mx' < 5 GeV for POMPYT, where Mx' includes the final 
state electron). We exclude these regions for the measurement of the diffractive structure func- 
tion by an upper cut on j3. To study acceptance and migration effects for these small masses 



we have generated an additional sample of exclusively produced p s [BH] 



The cuts given below to select diffractive events limit the acceptance for double-dissociative 



events, where the proton also dissociates. The PYTHIA Monte Carlo |33] has been used to 
study the detector response for the nucleon system M^ in double dissociation ■jp events. The 
nucleon system mass spectrum and the fraction of double-dissociative to single-dissociative 
events was taken from hadron data. 



Non-diffractive DIS processes were generated using the HERACLES 4.4 program 38 1 which 



incorporates first order electroweak corrections. The Monte Carlo generator LEPTO 6.1 |39 



interfaced to HERACLES via the program DJANGO 6.0 |^0|, was used to simulate QCD 
cascades and fragmentation. The parton cascade was modelled with the colour-dipole model 
including the boson-gluon fusion process by the ARIADNE 4.03 (CDMBGF) program. 



The fragmentation into hadrons was performed with the Lund string hadronisation model 
as implemented in JETSET 7.2 For the proton parton densities the MRSD'_ set K 



was 



chosen, which adequately represents our structure function results . In previous studies 

it was shown that this model gives a good description of the energy flow between the current 

jet and the remnant jet. 
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Combined sets of the Monte Carlo generators, described above, were used to simulate the 
expected final states in our DIS sample: the first one to describe non-diffractive DIS processes 
and the second one to model diffractive events. 

QED radiative processes were not simulated for diffractive events; however, with the selection 
cuts of section |, radiative corrections to the DIS cross sections are below 10% [0 and are 
expected to be of the same order for diffractive processes. All Monte Carlo events were passed 
through the standard ZEUS detector and trigger simulations and the event reconstruction 



package |jTO[- According to Monte Carlo studies, the efficiency of the trigger and of the final 
selection cuts is the same for diffractive and non-diffractive DIS events. The overall trigger 
acceptance is above 95%, independent of x and in the range of interest for this analysis. 



6 Event selection 



The selection of DIS events was similar to that described in our earlier publications |T8[ . 
The following offline cuts were applied: 

• Eg > 5 GeV, to ensure good electron identification; 

• QIa>8 GeV2; 

• yjB > 0.04, to give sufficient accuracy for DA reconstruction; 

• 5 > 35 GeV (with respect to the measured interaction point), to reduce radiative correc- 
tions and photoproduction background; 

• He < 0.95, to reduce photoproduction background; 

• the impact point of the electron on the face of the RCAL was required to lie outside a 
square of side 32 cm centred on the beam axis (box cut), to ensure that the electron 
shower was fully contained within the calorimater and its position could be reconstructed 
with sufficient accuracy; 

• a vertex, as reconstructed from VXD+CTD tracks, was required with \Zjjt^\ < 40 cm. 

In addition algorithms were used to reject cosmic-ray induced events and QED Compton events. 

A total of 31k events was selected in this way corresponding to an integrated luminosity of 
0.54 pb~^. Using the number of events produced by unpaired electron and proton bunches, 
the contamination from beam-gas background and from cosmic-ray muons were estimated to 
be less than 1% each. The background in the total DIS sample due to photoproduction was 
estimated to be (2.5 ± 1)% from a fit to the shape of the S distribution before the above cut 
on 6 was applied |TI 



7 Properties of diffractive events 

In the following section, the observed data distributions are compared to distributions from 
various Monte Carlo models. We discuss the criteria used to select the diffractive events, the 
methods used to determine their relative contribution and the observed inclusive distributions 
of the diffractive sample. 
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7.1 Selection criterion 



A large fraction of diffractive processes at HERA exhibit a rapidity gap in the main detector 
between the scattered proton system and the hadronic activity generated by the dissociation of 
the photon, while large rapidity gaps are suppressed in non-diffractive DIS events. Therefore 
the presence of a rapidity gap has been used as a selection criterion ||T], 0. An improved criterion 
to separate diffractive from non-diffractive events is presented here, which uses the direction of 
the total hadronic energy flow of the event, determined from all the detected particles in the 
final state. 

We define the maximum pseudorapidity of an event, r^max, as the maximum value of the pseu- 
dorapidity of all calorimeter condensates with energy greater than 400 MeV or tracks with 
momentum of at least 400 MeV/c. A condensate is a contiguous energy deposit above a mini- 
mum energy threshold. We studied the effect of varying the minimum energy -Emin = 400 MeV 
and found that above -Emin = 200 MeV the r^max distribution does not change significantly. We 
chose -Emin = 400 MeV as a conservative compromise between accepting diffractive events and 
rejecting noise. 

For values of ?7max up to 1-1.5 the non-diffractive DIS background is a negligible background 
to the diffractive sample, which increases for values of r/max above 1.5-2. In previous ZEUS 
publications [|^, 0] diffractive events were selected by ?7inax < 1-5. This cut selects a rather pure 
sample of diffractive events, useful to establish a signal but it limits acceptance for events with 
large Mx- 

The r^max-cut is dependent on the most forward condensate but does not use the information 
from the full energy flow. Larger acceptance can be achieved by including more information 
from the hadronic energy flow. Since in diffractive scattering the proton remains intact or, in 
the case of double-dissociative events, dissociates independently from the photon, the hadronic 
activity in the detector in general will not follow the proton direction. The hadronic angle 
Oh defined in eq. (|l|) represents the average direction of the hadronic activity. Non-diffractive 
DIS events have mostly cosOh near 1 because of the colour flow between the struck quark and 
the outgoing proton system, while a substantial fraction of diffractive events is found at cos9h 
less than 1. Figures b show scatter plots of ?7niax versus 9h for the diffractive and non- 
diffractive DIS Monte Carlo samples. A cut cos6h < 0.75 combined with r^max < 2.5 allows a 
larger acceptance of diffractive events than the ?7max < 1-5 cut, at the price of a slightly higher 
background which has to be subtracted. We call this combined cut, used to select the diffractive 
sample, the r/max-^'/f cut. 

Note that the term "diffractive" is used to indicate single diffractive dissociation of the photon 
together with that fraction of events where both the photon and the proton dissociate and the 
proton system is not detected. From proton-proton measurements of the ratio of double- to 
single-dissociative events, we estimate this ratio to be approximately 0.76 in the measured W 
range. As shown in Fig. |], we find that excited proton states with mass Mjy ^ 4 GeV would 
pass the diffractive selection cuts. Beyond this range the energy deposition in the forward 
calorimeter is typically above 400 MeV. The overall acceptance for double-dissociative events 
is 23%. We therefore have an estimated double-dissociative contribution of ~ (15±10)% which 
is expected to be independent of (3 and and not vary significantly with x^. This result 
assumes factorisation in Regge theory, i.e. that the nucleon mass spectrum and the ratio of 
double- to single- dissociative events is similar to that measured in proton-proton collisions at 
similar energies. 
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7.2 Estimation of the difFractive component 

In this section, only the shapes of the distributions and not the absolute normalisations of 
the diffractive models are considered. The ?7max and 6h distributions are used to determine 
the fraction of diffractive events passing the DIS selection criteria. A linear combination of 
diffractive DIS (NZ or POMPYT) and non- diffractive DIS Monte Carlo events are fitted to the 
data. 

The ?7max and Oh distributions were first fitted separately to check consistency between the 
results and then together to obtain a global result. Figure § shows the fits to these distributions. 
The part of each distribution that corresponds to the forward region of the detector (high values 
of ?7max and low values of Oh) was put into one single bin to reduce problems associated with a 
detailed description of the hadronisation of the proton remnant. For each distribution a variety 
of different binnings was tried and the results were found to be stable. 



Model 


^max 


Oh 


^max+^/f 




% of diflFr. 




% of diffr. 


Xdof 


% of diffr. 


Xdof 


NZ 


14.2 ± 2.5 


4.7 


15.3 ± 2.5 


4.0 


14.8 ± 3.0 


4.3 


SP 


35.9 ± 7.0 


10.5 


15.4 ± 2.6 


5.0 


33.0 ± 6.0 


24 


HP 


10.3 ± 2.0 


2.2 


10.8 ± 2.0 


4.1 


10.5 ± 2.7 


3.0 


HP+SP 


15.6 ± 1.3 


3.7 


13.4 ± 1.3 


5.4 


14.6 ± 1.4 


4.6 



Table 1: Fraction of diffractive events and per degree of freedom (Xdo/) values obtained from fits 
using NZ, SP, HP or HP+SP. 

The results are summarised in Table 1, where the default parameters have been used for the 
models. Since neither NZ nor POMPYT describes diffractive vector meson production a simu- 
lation of exclusive p° production was added in order to incorporate the effect of low-mass states. 
This contribution was estimated to be typically ~ 7% of the diffractive sample from a fit to the 
observed Mx spectrum in different intervals. For each model, a reduction of Xdo/ by 1-2 
was found when p° production was included, with consistent results obtained for the fraction 
of diffractive events. 

The SP model was also extensively tested. In fits to the r^max and Oh distributions, SP does 
not reproduce the shapes correctly. Its very soft [3 distribution tends to populate large r^max 
bins and, consequently, the fits do not describe the data. The inconsistency of the results 
obtained by the fits to the r^max and the Oh distributions shown in Table 1 indicate that a 
pure soft (3 distribution cannot describe the data. For these reasons the SP model alone is not 
considered any further. Results obtained with the combined HP+SP model, discussed in the 
following section, are also given in Table 1. The fractions obtained with the HP+SP model are 
similar to those determined using the NZ model. 

It is possible to explain the different predictions from NZ and HP models in terms of the 
j3 distribution used in these models: they both contain the hard component responsible for low 
??max (high Oh) events but the NZ model also contains a soft contribution which is predicted 
to be ~ 40% of the diffractive cross section. Most of the events originating from this soft 
component are hidden under the large background from normal DIS events and so are not 
accessible to the present study. 
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7.3 Inclusive distributions 



In the following, the shapes of the observed distributions in W, Q^, x, Mx, Xjp and (5 are 
considered. The relative normalisation of the models is obtained from the above fits. It should 
be noted that the normalisation of the non-diffractive component, which is relevant for the 
background subtraction, is independent of the diffractive model used to fit the data to within 
5%. 

In order to confine the analysis to regions of acceptance above ~ 80%, the following {Mx,y) 
intervals were considered: 

Mx < 10 GeV for 0.08 < y < 0.2 
Mx < 16 GeV for 0.2 < y < 0.3 
Mx < 20 GeV for 0.3 < y < 0.8 

According to Monte Carlo studies, the ?7max-6'iy cut reduces the non-diffractive DIS component 
by ~ 60% and the diffractive component by ~ 20%, giving a contamination from non-diffractive 
DIS of less than 15% in these {Mx, y) intervals. This background is subtracted from the data 
before comparison with the diffractive Monte Carlo predictions. 

Figure ^ shows the x, Q^, W, x^p, Mx and /? distributions after applying the rj^s^x-On cut, 
requiring the data to be in the accepted ranges of (Mx, y) and subtracting the DIS back- 
ground indicated in the figure. The errors on the data points are calculated by summing in 
quadrature the statistical error (which is the dominant error) and 50% of the total subtracted 
DIS background (which is taken as a conservative estimate of the uncertainty due to the DIS 
background). In addition, the predictions from the two diffractive models (NZ and HP) are 
shown. 

In general, both models describe the data. Differences are observed in the Mx and (3 distribu- 
tions, where the HP model underestimates the observed number of events at low (3 values and 
does not reproduce the observed Mx distribution at large Mx- The NZ model, incorporating 
a soft component, describes the observed /? and Mx distributions. 

A pure "hard" (3 distribution cannot account for the data, therefore, the observed /3 spectrum 
was fitted as a sum of a "hard + soft" contribution from the POMPYT Monte Carlo. This 
resulted in a contribution of ^ 60% and ~ 40% from HP and SP, respectively. This HP+SP 
model is also shown in Fig. ^ Comparison with the data indicates that such a model also 
describes the observed jS behaviour. 

To investigate the f3 distribution in more detail, each {Mx,y) interval was divided into two Q'^ 
bins: 

g2 = 8 _ 20, 20 - 160 GeV^ 

The results, together with the predictions from the diffractive models are shown in Fig. ^ In 
general, the Monte Carlo models reproduce the shape of the data reasonably well. However, 
in the high-y and low-Q^ intervals where the mass extends to larger values (Fig. ^a and c), 
the soft contribution is important. The NZ model describes the data best in this region. The 
HP+SP model reasonably describes the data and gives an improved description compared to 
the HP model in each {Mx,y) interval. 
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Using the NZ model, the combined fit to the r^max and 6h distributions was performed in bins of 
W and X respectively, separately for the two Q'^ intervals indicated above, to extract the fraction 
of diffractive events as a function of these variables. Figure |^ shows the diffractive fraction as a 
function of W and x for different values of Q^. The results extracted using the HP+SP model 
agree within statistical errors. The results extracted using the HP model give a normalisation 
which is ~ 30% lower, but with the same dependence on x, W and Q^. The fits are mainly 
sensitive to the hard component: a large uncertainty on the diffractive contribution to the DIS 
sample comes from the soft part in the pomeron structure function, which is suppressed by the 
applied cuts, especially at small values of W. In all cases, no strong dependence of this ratio is 
observed as a function of x, W or Q"^. 



8 Measurement of the diffractive structure function 

As described in section 0, the differential cross section can be expressed in terms of the diffractive 
structure function as a function of (3, x^ and Q^. In this section, we discuss the resolution 
of the measured quantities and describe the kinematic region chosen. We finally discuss the 
systematic errors and present the results of the measurement of F^^^\ 

8.1 Extraction of ^2^^^^ 

According to Monte Carlo studies (see section |^), the resolution of is 25%, independent 
of Q^. The resolution of x varies smoothly with x from 20% at x = 10"^ to 50% at 10~^, 
almost independent of Q^. The resolution of Mx, reconstructed with the method described in 
section |, is approximately 27%, independent of Mx- The Mx reconstruction is affected by 
energy loss in inactive material in front of the calorimeter and the position determination of 
hadrons. In order to reduce migrations at small masses, the cell energy thresholds for isolated 
cells were increased. Monte Carlo studies show that, except for very small masses (< 3 GeV) 
where calorimeter noise becomes important, Mx is systematically shifted by 10% to smaller 
values, independent of y and Q^. In order to compensate for this shift, a correction factor of 
1.10 was applied to the measured Mx values for the determination of the diffractive structure 
function. The resolution of Xp is approximately 25%. The resolution of (3 varies smoothly with 
(3 from 40% at (3 = 0.1 to 20% at /3 = 0.8. 

Below of 8 GeV^, the event acceptance drops below 50% due to the box cut requirement. 
The statistics of the 1993 data allow four ranges in to be selected above this lower limit. 
The migration of events is large at small values of x: we therefore chose bins where the central 
x-value is above 4 ■ 10~^. The acceptance of the diffractive component increases as a function 
of y: we therefore select only bins with y > 0.08. The overall acceptance due to the DIS and 
diffractive cuts in the selected bins given in Table 2 is always greater than 50% and typically 
~ 80%. The Mx resolution determines the chosen bin size in the variables f3 and x^. The 
purity, defined as the fraction of simulated events generated in a bin and measured in the same 
bin, is always greater than 25% and typically ~ 40% in each of the selected bins. 

In order to control the influence of photoproduction background, radiative corrections and Fi 
contributions, we restrict our analysis to y < 0.5. As a consequence the minimum scattered 
electron energy requirement is raised to 10 GeV. We checked that our sensitivity to Fl is smaller 
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than the quoted errors in all bins. Furthermore, the region (3 < 0.8 is selected to exclude the 
region of low masses where vector meson production is dominant. 

The level of photoproduction background is estimated in bins of Xp and by fits to the 5 
distributions (see [1^ for details). Since it is typically ^ 1% and always below 4% we do not 



correct for this background. 

We select bins with < 0.01 and (3 > 0.1 where the non- diffract ive component can be safely 
estimated. In each of the bins the number of events is then evaluated by subtraction of the 
estimated number of DIS background events, based on the ARIADNE Monte Carlo program 
with the normalisation described in section 6. The contribution of the DIS background is given 
in Table 2. 

To unfold the effects of acceptance and event migration we used the NZ Monte Carlo event 
sample, which gives a good description of our data. For this initial study we used a one-step 
matrix unfolding procedure and applied a bin-centring correction for the quoted F^^^^ values. 



8.2 Systematic errors 



Several systematic checks were performed to estimate the uncertainties due to the selection 
cuts, background estimate and the unfolding. Systematic errors due to the DIS event selection 
were evaluated in the following way (see [|TB] for a detailed discussion): 



• different algorithms were used to identify the scattered electron which differ in purity and 
efficiency. The changes to ^2^^^'' were below 10%; 

• the cut on E'^ was decreased from 10 to 5 GeV to study the effect of a possible mismatch 
of the shower profiles of data and Monte Carlo at small energies. The change of F^^^^ 
was less than 5% in each bin; 

• the box-cut was changed by 2 cm from the nominal values, to study the effects of electron 
position reconstruction at small angles. This resulted in changes which were always less 
than 15%; 

• the 5-cut was raised from 35 GeV to 40 GeV, to study the effect of radiative corrections, 
which were not included in the simulations. This resulted in a general shift of ^ 10% 
towards smaller F^^^^ values; 

• the yjB-cnt was changed from 0.04 to 0.02 and to 0.06. This affected the region of large 
Xjp where F^^^^ changes by about 10%. 



Systematic errors due to the diffractive event selection were evaluated in the following way: 



• the effect of a possible mismatch between the hadronic energy scale in the Monte Carlo 
and the data was investigated by shifting the hadron energy scale by 7% in the Monte 
Carlo simulation. The use of the DA variables resulted in changes on F^^^^ which were 
always smaller than 2%; 
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• the fraction of low-mass events was reduced by 50%. Due to migrations from (3 > 0.8, 
this change influences the small Q^, high f3 bin, where the values were shifted upwards 
by ~ 10%; 

• the HP model was used instead of the NZ simulation for unfolding the data. Some effect 
was seen in the small /9-region, where the pomeron structure functions differ. The changes 
to ^2^''^^ were typically ~ 10%; 

• as a systematic check for the estimate of the DIS background the ?7max-cut was reduced 
from 2.5 to 2.0 resulting in changes of up to 20% in the highest bins. The //max-cut 
was also increased from 2.5 to 3.0 resulting in changes of up to 10%; 

• similarly, the 6h cut was removed, yielding changes below 5%; 

• the cells with rj > 2.5 were removed to check the dependence on the double-dissociative 
contribution, resulting in changes which were up to 5%. 

Overall most of these checks yielded results which agree with the standard method within 
statistical errors. The differences of the DIS and diffractive systematic checks compared to the 
standard method were combined in quadrature to yield the quoted systematic errors. 



8.3 Results 

Table 2 summarises the results for F^^^\ for the 0.54 pb~^ (±3.5%) integrated luminosity. The 
statistical errors include statistical uncertainties from the Monte Carlo models used for the 
unfolding. 

The ^2^*'^'* results are displayed in Fig. The data are observed to fall rapidly as a function of 
increasing x^. In the measured bins the dependence of F^^^^ on at fixed (3 values is weak. 
We have investigated whether the x^-dependence of F-^'^^^ is the same in each (3, interval, 
as expected if factorisation holds. For this purpose we performed fits of the form: 

bi ■ {i/x^r 

where the normalisation constants bi were allowed to differ, while the exponent was the same 
for each (3, interval. The result of the fit was: 

a = 1.30 ± 0.08 (stat) + (sys). 

The systematic errors are calculated by re-fitting the ^2^''^^ values according to the variations 
listed in section and combining the positive or negative deviations from the central value 
of a in quadrature. The overall statistical values of these fits are in the range 8.2-14.0 
for 23 degrees of freedom depending on the systematic check. The values for each of the 
f3,Q^ intervals are in the range 0.1-1.1 per degree of freedom. Within the present accuracy, 
the data are therefore consistent with the assumption of factorisation in the measured kine- 
matic range. The value of a is consistent with recent results from the HI Collaboration of 
a = 1.19 ±0.06 ±0.07 ||. 
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The observed dependence on Xp is steeper but still compatible with a Donnachie-Landshoff 
type of flux factor which yields a ^ 1.09 and which is based on a phenomenological de- 
scription of "soft hadronic" diffractive interactions. The modified flux of Goulianos yields 
a {l/XjpY dependence with a ~ 0.93, a value which is disfavoured by the data. 

In order to illustrate the (3 and dependence of F^^^\(3, Q"^, x^), we integrated F^^^^ over the 
measured range of x^, 6.3 ■ 10~^ < x^ < 10~^, using the fitted x^ dependence. The resulting 
values of F^{P, Q"^) are shown in Fig. ^as a function of [3 and Q^. It should be noted that these 
results assume that a universal x^ dependence holds in all regions of (3 and . In particular, 
there is a contribution due to regions of x^ which are not measured and where the hypothesis 
of a universal Xjp dependence has not been tested experimentally. 

The F^{P,Q'^) values as a function of P for fixed Q"^ are consistent with a fiat /? dependence 
as expected from the aligned jet model |^ and the model of Buchmiiller As a function 



of for fixed /3, the F<P{P,Q^) values are approximately independent of for all /? values, 
which is consistent with a picture where the underlying interaction is the scattering of a virtual 
photon with point-like quarks within the pomeron. 

As a next step we determined a compact parametrisation for the F-^'^^^ results, which is also 
shown in Fig. |] and Fig. where the following form was adopted: 

F,''^'^ = {i/x,r-b.m-f3) + ^-{i-(3n 

with a = 1.30. This parametrisation assumes factorisation and no Q^-dependence. The soft 
contribution to the structure function was considered by the inclusion of the (1 — term. 
The multiplicative factor of | was chosen such that the integral over (3 of the soft contribution 
is equal to that of the hard contribution when c = 1. The power of 2 was adopted from the NZ 
model; this assumption cannot be tested with the current measurement. The results of the fit 
were: 

b = 0.018 ± 0.001 (stat) ± 0.005 (sys), 

c = 0.57 ± 0.12 (stat) ± 0.22 (sys), 

with a statistical the range 15-23 for 33 degrees of freedom depending on the systematic 
check. A fit without the (1 — soft contribution resulted in values in the range 56-81 for 
34 degrees of freedom. This increased value indicates that a soft component is required in 
the pomeron structure function. 



We now consider the cross section predictions of the models discussed in section ^. The Fr 



D(3) 



results are displayed in Fig. [1^ where the data are compared with the predictions of several 
models of single-diffractive dissociation for which the momentum sum rule for quarks is not 
satisfied. The estimated 15% fraction of double-dissociative events has been subtracted in order 
to compare with these models. 

At high /3-values the predictions of Nikolaev-Zakharov, Donnachie-Landshoff and Capella et 
al. underestimate the observed values slightly, but are generally in reasonable agreement. At 
smaller /3-values, the Donnachie-Landshoff parametrisation, which includes only a hard com- 
ponent of the pomeron structure function, underestimates the observed F^^^\ The Capella 
et al. and Nikolaev-Zakharov predictions, which also include a soft component, are able to 
give a fair description at smaller /5-values. The factorisation-breaking effects in the model 
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of Nikolaev-Zakharov, which occur at small (3 values, are too small to be observable in this 
analysis. 

In Fig. |lT] the data are compared with a model for which the momentum sum rule for the 
pomeron structure function is assumed for the light quark flavours (u,d) and the (3 dependence 
is taken from the parametrisation, discussed above. Adopting the Donnachie-Landshoff flux 
factor, the observed F^^^^"* is about a factor three to four below the expectation if the momentum 
sum rule is assumed to be fulfilled only by quarks. An uncertainty arises from the choice of 
the pomeron flux factor: if the Ingelman-Schlein form for the flux factor is adopted then the 
prediction is reduced by approximately 30%. Even if the Goulianos prescription for the flux 
is adopted, the observed ^2^^^"* results are always below the predictions. These comparisons 
indicate that in an Ingelman-Schlein type model the quarks alone inside the pomeron do not 
satisfy the momentum sum rule. 



9 Conclusions 

The properties of diffractive DIS events with > 8 GeV^, selected by a large rapidity gap 
requirement, have been investigated. Different Monte Carlo models, such as the POMPYT 
model, with a soft plus a hard pomeron structure function, or the Nikolaev-Zakharov model 
describe the shape of the observed kinematic distributions. Using these models, reliable ac- 
ceptance corrections for the measured data can be obtained and corrected cross sections can 
be determined. The relative contribution of diffractive events to the total DIS cross section is 
found to have no strong dependence on x, W oi . 

The diffractive proton structure function ^2^*^^'' is presented, integrated over t, the square of the 
momentum transfer at the proton vertex, as a function of x^, the momentum fraction lost by 
the proton, of (3, the momentum fraction of the struck quark with respect to x^, and of Q^. The 
structure function is measured in the kinematic range of 0.08 < y < 0.5, 8 < < 100 GeV^, 
6.3 • 10~^ < Xjp < 10~^ and 0.1 < /5 < 0.8. Within the experimental errors, the measurement is 
consistent with models where diffraction is described by the exchange of a particle-like pomeron 
where the structure function factorises into a pomeron flux factor, which depends on and 
a pomeron structure function, which is independent of x^p. The diffractive structure function 
is also well-described by the Nikolaev-Zakharov model, which does not require the concept of 
a particle-like pomeron, in terms of overall normalisation and dependence on the kinematic 
variables, x^, [3 and Q^. The x^ dependence is consistent with the form {1/XjpY where 
a = 1.30 ± 0.08 (stat) (sys) in all bins of /3 and Q^. The value of a is slightly higher but 

compatible with that obtained from hadron-hadron interactions and in agreement with recent 
results from the HI collaboration. In the measured range, the pomeron structure function 
is approximately independent of at fixed (3 consistent with an underlying interaction where 
the virtual photon scatters off point-like quarks within the pomeron. The /3-dependence of the 
pomeron structure function requires both a hard and a soft component. In an Ingelman-Schlein 
type model, where commonly used pomeron flux factor normalisations are assumed, it is found 
that the quarks within the pomeron do not saturate the momentum sum rule. 
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Table 2: ZEUS 1993 F2 results. The overall normalisation uncertainty of 3.5% is not in- 
cluded. The data contain an estimated 15±10% fraction of double-dissociative events. 
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Figure 1: Diagram of a difFr active event. 
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Figure 2: versus rymax distribution for diffractive (NZ) and non-diffractive (CDMBGF) Monte 
Carlo events and for the selected DIS data. The full line indicates the 77niax-^if cut used to select 
diffractive events. The dotted line corresponds to ?7niax=l-5. 
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Figure 3: Acceptance for double dissociative events. The mass of the nucleon system, Mn, for double 

dissociative events generated by the PYTHIA Monte Carlo is indicated by the full line histogram. 
The shaded area indicates those events which are selected by the ??max-^H cut. The fraction of double 
dissociative events rejected by this cut, as a function of Mn, is indicated by the line. 
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Figure 4: Percentage of DIS data as a function of r/max and 9h- The data are described by the sum of 
the diffractive and non-diffractive contributions obtained from Monte Carlo simulation, with relative 
fractions determined by a fit to the data. The dashed line corresponds to the diffractive contribution 
and the sum of the diffractive and non-diffractive Monte Carlo models is indicated by the full line. 
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Figure 5: Observed distributions of x, Q"^, W, x^, Mx and P for the selected {Mx,y) intervals. 
Uncorrected data are indicated by the dots. The errors are the statistical errors combined in quadrature 
with 50% of the non-diffractive DIS background. The predictions from HP (full line), HP+SP (dashed 
line) and NZ (dotted line) models are shown. The non-diffractive DIS background which has been 
subtracted from the data is indicated by the shaded area. 
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Figure 6: Observed /3 distribution as a function of (y, Q'^,Mx)- The intervals are 8-20 and 20-160 
GeV^, the y intervals are 0.08-0.2, 0.2-0.3 and 0.3-0.8, and the Mx intervals are (a,b) 0-20, (c,d) 0-16 
and (e,f) 0-10 GeV. Uncorrected data are indicated by the dots. The errors are the statistical errors 
combined in quadrature with 50% of the non-difFractive DIS background. The predictions from HP 
(full line), HP-I-SP (dashed line) and NZ (dotted line) models are shown. The non-diffractive DIS 
background which has been subtracted from the data is indicated by the shaded area. 
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Figure 7: Observed fraetion of diffractive events as a function of Wda and xjja in two intervals. 
The data are fitted to the NZ model for diffractive processes and the CDMBGF model for the non- 
diffractive contribution. The errors are the statistical errors combined in quadrature with 50% of the 
non-difFractive DIS background. 
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Figure 8: The results of F2 ^^'{l3,Q^,Xjp) compared to the parametrisation discussed in the text. 
The inner error bars show the statistical errors, the outer bars correspond to the statistical and DIS 
event selection systematic errors added in quadrature, and the full line corresponds to the statistical 
and total systematic errors added in quadrature. Note that the data include an estimated 15% 
contribution due to double dissociation. The overall normalisation uncertainty of 3.5% due to the 
luminosity uncertainty is not included. 
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Figure 9: The results of F2^{P, Q^) compared to the parametrisation discussed in the text, indicated 
by the full line, and the — /3) hard contribution, indicated by the dashed line. The inner error bars 
show the statistical errors, the outer bars correspond to the statistical and systematic errors added in 
quadrature. The systematic errors combine in quadrature the fits of the Xjp dependence due to each 
of the systematic checks discussed in the text. Note that the overall normalisation is arbitrary and 
is determined by the experimental integration limits over (6.3 • lO"-^ < Xjp < 10-2). The data 
include an estimated 15% contribution due to double dissociation. 
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Figure 10: The results of -F2 compared to various models discussed in the text. Note that the 
estimated 15% contribution due to double dissociation has been subtracted in order to compare with 
models for the single dissociation cross section. The inner error bars show the statistical errors, the 
outer bars correspond to the statistical and DIS event selection systematic errors added in quadrature, 
and the full line corresponds to the statistical and total systematic errors added in quadrature. The 
overall normalisation uncertainty of 3.5% due to the luminosity and 10% due to the subtraction of the 
double dissociation background is not included. 
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Figure 11: The results of F2 compared to an Ingelman-Schlein type model for which the momentum 
sum rule (MSR) for quarks within the pomeron is assumed. The f3 dependence is taken from the 
parametrisation discussed in the text. Note that the estimated 15% contribution due to double 
dissociation has been subtracted in order to compare with models for the single dissociation cross 
section. The inner error bars show the statistical errors, the outer bars correspond to the statistical 
and DIS event selection systematic errors added in quadrature, and the full line corresponds to the 
statistical and total systematic errors added in quadrature. The overall normalisation uncertainty of 
3.5% due to the luminosity and 10% due to the subtraction of the double dissociation background is 
not included. 
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